function A = assemble_matrix_1D(?)


A = sparse(matrix_size(1),matrix_size(2));

for n=1:number_of_elements
   
    for alpha=1:number_of_local_basis_fun_trial
        
       for beta=1:number_of_local_basis_fun_test
          
           int_value= Gauss_quad_1D_trial_test(?);
           
%          i = Tb(beta,n);
%          j = Tb(alpha,n);
           T(Tb(beta,n),Tb(alpha,n))=T(Tb(beta,n),Tb(alpha,n))+int_value;
           
       end
        
        
    end
    
    
    
    
    
    
end